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Abstract 



We conjecture that a new class of 3d Af = 2 theories are associated with a quiver Q 
and a mutation sequence m on it. We define the cluster partition function from the pair 
(Q, m), and this partition function coincides with the S$ partition function of the associated 
3d Af = 2 theory T[(Q, m)]. 

Our formalism includes the case where 3d Af = 2 theories arise from the compactification 
of the 6d (2,0) Ajy_i theory on a large class of 3-manifolds M, including complements of 
arbitrary links in S 3 . In this case the quiver is defined from a 2d ideal triangulation, the 
mutation sequence represents an element of the mapping class group, and the 3-manifold is 
equipped with a canonical ideal triangulation. Our partition function coincides with that of 
the holomorphic part of the SL(N) Chern-Simons partition function on M; when N = 2 and 
M is hyperbolic, the partition function reproduces the gluing conditions of ideal hyperbolic 
tetrahedra in the semiclassical limit. 
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1 Introduction 



It has recently been discovered [1, 2, 3] (see also earlier works [4, 5, 6]) that there exists a 
beautiful correspondence ("3d/3d correspondence") between the physics of 3d M = 2 gauge 
theories and the geometry of (hyperbolic) 3-manifolds. The latter, in more physical language, 
is the study of analytic continuation of 3d pure Chern-Simons SU(2) gauge theory into a 
non-compact gauge group SL{2) [7, 8]. More quantitatively, one of the consequences of this 
correspondence is that given a 3-manifold M there is a corresponding 3d M = 2 theory T[M] 
such that the partition functions of the two theories coincide [1]: 

Z*- d { $f 2 [S$] = Z Chcrn - Simons [M] , (1.1) 

where the left hand side is the 3-sphere partition function [9, 10, 11] with a 1-parameter 
deformation by b [12], and the right hand side is the holomorphic part of the SL(2) Chern- 
Simons theory with level t. The two parameters b and t are related by b 2 ~ 1/t. 1 2 

The correspondence (1.1) provides a fresh perspective on the systematic study of a large 
class of 3d M = 2 theories, and relations between them. An arbitrary hyperbolic 3-manifold 
could be constructed by gluing ideal tetrahedra, and correspondingly we could construct 
complicated 3d Af = 2 theory starting from free M = 2 hypermultiplets. Gluing ideal tetra- 
hedra is translated into the gauging of the global symmetries, and the change of polarization 
into a the Sp(2N, Z) action [20, 21] on 3d M = 2 theories. The 2-3 Pachner move, which 
represents the change of the ideal triangulation of the 3-manifold, is translated into the 3d 
M = 2 mirror symmetry [2] . 

Despite the beauty of this correspondence, we should keep in mind limitations of the 
correspondence (1.1) - - not all 3d Af = 2 theories are of the form 7~[M]. The natural 
question is whether we can generalize the correspondence to a larger class of 3d M = 2 
theories beyond those associated with 3-manifolds, or more generally whether there are any 
geometric structures in the "landscape" or "theory space" of 3d M = 2 theories. To answer 
these questions it is crucial to extract the essential ingredients from the correspondence (1.1). 

Our answer to this question is that it is the mathematical structures of quiver mutations 
and cluster algebras which is essential for the correspondence (1.1). We study 3d M = 2 
gauge theories T[M] for a large class of hyperbolic 3-manifolds, including arbitrary link 
complements in S* 3 . 3 The present paper generalizes the previous works on this subject by 

1 Note that in (1.1) the same data appears in rather different forms on the two sides. For example, the 
geometry of M for the right hand side determines the choice of the theory T[M] itself on the left. Similarly, 
the deformation of the geometry, the parameter b, on the left hand side is translated into a parameter of the 
Lagrangian on the right. 

2 The first evidence for this conjecture [1] came from a chain of arguments involving quantum Liouville and 
Teichmuller theories. The semiclassical (t — > oo) expansion of the right hand side (1.1) reproduce hyperbolic 
volumes and Reidemeister torsions of 3-manifolds [13, 14]. See [15, 16, 17, 18, 19] for further developments 
in the 3d/3d correspondence. 

3 Suppose that we have a 3-manifold M, and a link L inside. The complement of a link is defined as 
a complement of a thickened link. More formally, the link complement is a complement of the tubular 
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the authors [1, 13], and is a companion to the previous paper with K. Nagao [14]. It is also 
closely related with [3] (see also [22, 19]) 4 . 

Our approach turns out to be much more general than (1.1), and includes theories asso- 
ciated with SL(N) Chern-Simons theories on 3-manifolds (cf. [23]), or more general theories 
not associated with 3-manifolds (Figure l). 5 We define a 3d M = 2 theory T[(Q,m)] for a 
pair of a quiver Q and a mutation sequence m on it, 6 satisfying the relation 

y3d Af=2 rc3l ^cluster /i r>\ 

Z T[(Q,m)]PfcJ - Z (Q,m) > 

where the right hand side is the cluster partition function defined in this paper. The right 
hand side contains a quantum parameter q, which is related to the parameter b on the left 
hand side by the relation (2.6). It will an exciting possibility to explore the properties of 
these theories further, part of which will be the subject of a separate publication [24]. 

3d M = 2 theories 

T[(Q,m)] 

Figure 1: The class of 3d M = 2 theories (denoted by T[(Q, m)}) is a special subset of general 
3d M = 2 theories, and is large enough to contain theories T[M] dual to the geometry 3- 
manifolds. Little is known at present about the detailed properties of the theories T[(Q, m)], 
when the theory cannot be written in the form T[M]. 

Basic Idea 

Let us here briefly summarize the basic idea behind our algorithm of identifying a 3d M = 2 
theory from a given 3-manifold, following closely the approach of [1]. 

neighborhood N(L) of M, i.e, M\N(L). By construction the boundary of the link complement is a disjoint 
union of 2d tori. 

4 However it is worthing emphasizing that their braid (branched locus) is not our braid, see further 
comments in section 3.1. 

5 All our theories are contained in the theories of "class 1Z" in [16]. However, for more general theories 
not described by quiver mutations there is no simple answer to For comparison one might be tempted to 
call the theories T[M] to be of "class Ai" (Ai for manifold) and theories T[(Q,m)] to be of "class C" (C for 
cluster algebras). 

6 The theory in addition depends on the choice of the boundary condition, as will be explained in section 
4.2. 
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Let us first consider a 3-manifold of the form Ex/, where £ := S ffi ^ is a Riemann 
surface with ft, punctures and genus g, and I is an interval of finite size (Figure 2). Such 
a 3-manifold is called a mapping cylinder. We will later generalize our analysis to more 
general 3-manifolds. We assume E is hyperbolic, i.e., x(^) < 0- When the surface £ has 
punctures, the trajectory of the punctures sweeps out the Id defect inside the 3-manifold, 
defining braids inside M (Figure 3). 




Figure 3: The motion of the punctures represents braids inside the mapping cylinder. 



We consider SL(2) Chern-Simons theory 7 on this 3-manifold (Figure 2). 

Since there is a canonical preferred direction on this 3-manifold, we could regard that 
direction as a time direction and carry out canonical quantization, obtaining the Hilbert 
space of our theory %y. on our Riemann surface £. This is identified with the Hilbert space 
of the quantum Teichmiiller theory, formulated first in [25, 26]. 8 

7 We discuss the holomorphic part of the SL(2) Chern-Simons theory, and for this purpose it does not 
matter whether the gauge group is SX(2,C) or SX(2,R). 

8 We do not give a detailed explanation for this reason here, see the review material in [1]. Here it suffices 
to point out that the classical saddle points of 3d SX(2,R) Chern-Simons theory are given by flat SX(2,R) 
connections on the 3-manifold, whereas Teichmiiller space a connected component of moduli space of flat 
SX(2,R) (or rather PSL(2,M)) connections on E, and hence the latter is a naturally arises in the canonical 
quantization of the latter. 
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Since Chern-Simons theory is topological, the time evolution is trivial, and the only non- 
trivial information in this case is the choice of the boundary conditions. In other words, 
corresponding to the two boundaries we need to specify the two states 

| in), |out) G H s , (1.3) 

and the partition function is simply given by the overlap of the two states 

Z Sx/ = (in|out) . (1.4) 

The coordinate description of "Hs depends on a choice of an ideal triangulation, and the in 
and out states in (1.4) might be naturally described in different triangulations. To take this 
into consideration we introduce an operator representing the change of the triangulation 
between in and out states, leading to 9 

Z Ex/ = (in | | out) . (1.5) 

Geometrically ip will be an element of the mapping class group of E, i.e. tp is a large coordinate 
transformation on E. 

We can also identify the in and out states, and take a sum over all the possible states. 

%x5!) v = Tr(v?) = J d(ixi) (in\0 |in) . (1.6) 

In this case the geometry is that of a mapping torus (E x S 1 )^, defined by identifying 
(x, 0) ~ 1) for E x I. Here we have taken I = [0, 1], and ip to be an element of the 

mapping class group of E. 

Since the action of is given explicitly, we can evaluate this expectation value and 
obtain an integral expression — the integrand contains one quantum dilogarithm function 
[27, 28, 29] for each flip of the triangulation. We can then read off the corresponding 3d 
M = 2 gauge theory, using the relation (1.1) as a guideline. The resulting 3d M = 2 theory 
can be thought of as a theory on the duality domain wall theory [30, 31] inside the 4d M = 2 
theory of [32]. 

9 We can describe this more formally. Given two triangulations T, T", there is an operator if>T',T such 
that 

|out)T' = <Pt',t |out)r . 

and the partition function is given by 

Z^xi = (in| T |out) T / = (in| T (<Pt',t) |out) T . 

The difference between | )t and | )t> will be crucial when we identify the boundary data, see (1.6) and section 
3.2. 
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The Results 



While the strategy outlined to this point should work, this program has never been worked 
out in generality. It is also the case that the resulting 3-manifold is apparently limited to 
mapping cylinders or mapping tori, and it is not clear if this method generalizes to more 
general 3-manifolds. 

The goal of this paper is to fill in these gaps. In particular this paper contains the 
following results: 

1. The results of this paper are described nicely in terms of quiver mutations and cluster 
algebras (section 2). We introduce a new object, a mutation network (section 2.4), 
which encodes the combinatorial data of a quiver Q and a mutation sequence m. 

2. We will obtain an explicit integral expression for the partition function Z^q s ^ (cluster 
partition function) of (1.5), determined from the combinatorial data of the mutation 
network (section 2.3, in particular (2.30)). 

3. When the pair (Q,m) satisfies certain conditions, we can construct an associated 
mapping cylinder (E,</>) with a canonical ideal triangulation determined from (Q,m) 
(section 3.1). 

4. By appropriately identifying boundaries of the mapping cylinder (£,</?), we obtain a 
large class of hyperbolic link complements, including arbitrary link complements in S 3 
(section 3.2). This operation has a counterpart in the mutation network as well as the 
partition function. 

5. We construct the 3d M = 2 theory 1~[(Q, m)] whose S% partition function coincides 
with the partition function Z^q s ^ defined previously from the mutation network (1.2) 
(section 4.2). 

6. We consider dimensional reduction of our 3d M = 2 theory to 2d M = (2, 2) theory. The 
twisted superpotential coincides with the Neumann- Zagier potential [33] of hyperbolic 
geometry, and the vacuum equations of the 2d M = (2,2) theory reproduce the gluing 
equations of hyperbolic tetrahedra. 

The rest of this paper is organized as follows. In section 2 we derivate an integral 
expression of our partition functions, based on the formalism of quiver mutations and cluster 
algebras. The combinatorial data is summarized in the mutation network. We also write 
down the associated cluster partition function. In section 3 we explain the geometry of 
hyperbolic 3-manifolds and their ideal triangulations, associated with a mutation network 
satisfying certain conditions. In section 4 we study the 3d M = 2 theory dual to the 
mutation network. The final section (section 5) summarizes the results and comments on 
open problems. We include three appendices. Appendix A summarizes the properties of the 
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quantum dilogarithm function used in the main text. Appendix B summarizes hyperbolic 
ideal triangulations and gluing equations. Appendix C explains the effect of the Dehn twist 
on the triangulation. 

We have tried to make the paper accessible to a wide spectrum of readers, including 
mathematicians. In fact, most of the materials in sections 2 and 3 (apart from examples in 
section 3.4) require little prior knowledge of the subject (in physics and in mathematics), 
and no knowledge of supersymmetric gauge theories are necessary until section 4. Readers 
not interested in 3-manifold cases can skip section 3 and proceed directly to section 4. 

2 Quivers and Clusters 

In this section we summarize the construction of the Hilbert space %q associated with a 
quiver Q, and the action of the mutations m on "Hq. We also define the cluster partition 
function Z^q s ^, whose integral expression will be given in this section. 

This section will be formulated in terms of quiver mutations and cluster algebras [34] 
(see [35] for introduction). 10 While this might look abstract to some readers, this turns out 
to be the crucial mathematical structure underlying our 3d M = 2 theories. 

2.1 Quiver Mutations and Cluster Algebras 

Let us begin with a quiver Q, i.e., a finite oriented graph. We denote the set of the vertices 
of the quiver by /, and its elements by G I. 

In this paper, we always assume that a quiver has no loops and oriented 2-cycles (see 
Figure 4). This will always be satisfied for our quivers coming from ideal triangulations of 
E, discussed in next subsection. 



Figure 4: A loop (left) and an oriented 2-cycle (right) of a quiver. 

For vertices i,j G /, we define 11 

Qij := (jjarrows from % to j} — fj{arrows from i to j} , (2-1) 

i.e., \Qij\ represents the number of arrows from the vertices i to j, and the sign represents 
the chirality (orientation) of the arrow. Note that the quiver Q is uniquely determined by 
the matrix Qjj under the assumptions above. 

10 For appearance of cluster algebras in 4d Af = 2 theories, see for example [36, 22]. 
11 Qij here is denoted by Q(i,j) in [14]. 
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Given a vertex k, we define a new quiver fi k Q (mutation of Q at vertex k) by 



-Qi 7 (i = k or j = k) , , 

((J>kQ)ij ■= { (2.2) 

where we used the notation 

[x] + :=max(x,0) . (2.3) 

In more physical language, this can be regarded as a somewhat abstract version of the Seiberg 
duality [37]. However, the difference from there is that our mutations in general are outside 
the conformal window. 

Let us now construct a non- commutative algebra Aq associated with the quiver Q. This 
is generated by a set of variables y^, 12 associated with each edge i, satisfying the relation 13 

:= {yi (i6 j) | YjYi = <f Q Y.y,} , (2.5) 

where q is the quantization parameter which is related to the parameter b by 

q = e^ nb2 . (2.6) 

This parameter b will later be identified with the deformation parameter of S 3 in (1.1). The 
semiclassical limit is given by q — > 1, or b — > 0. The variables are the quantum versions 
[38, 39] of the so-called y- variables (coefficients) in the cluster algebra literature. 
We can promote the mutation \i k to an operator fi k on Aq: 

fa: A Q ^ A^q . (2.7) 



12 In section 4.2 we will comment on the interpretation of these variables as loop operators in 3d M = 2 
theories. 

13 More formally, we can think of this space as generated by y Q with a € Z^l, satisfying the relations 

g<«,/J) y «y - y Q+/J , (a, (3) = -(13, a) ^ f aQ(3. (2.4) 
In this notation yi = y £i for the i-th basis of Z' 7 ' . 
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The operator jl k acts on yj by 



y'i ■= Afe(yi 



yr 1 , (i = k) 

\Qm\ 



-,Qik[Qki]+^. Kll [Qki[+ ^_|_g-sgn(Q w )(2m-l)y^-Bgn(Qfci) 



m=l 



q Qik[-Qki]+y.y k l-Qki]+ JJ ^ + ? sgn(Q fci )(2m-l) Vfc -l^-sgn(Q fc i) 



m=l 



(2.8) 



We can naturally extend the action of ^ to the whole of Aq. The resulting variables y£ 
satisfy commutation relation (2.5) for fikQ, hence fik is indeed a map from Aq to A^q. 

The commutation relations (2.5), if written in terms of variables Y, = log(yj), take a 
simple form 

[Y ] ,Y t } = 27rb 2 V^lQ hj . (2.9) 

This has a standard representation on a Hilbert space — we can choose a polarization, i.e., 
perform linear transformations to find coordinate and momentum variables, and coordinates 
(momenta) act by multiplication (differentiation). 14 We denote this Hilbert space Hq; the 
algebra Aq is now the set of operators acting on this state. We will present more concrete 
discussion in section 2.3. 

In the following we consider a sequence of quiver mutations (/x mi , . . . , l^m L ), specified by 
m = (mi, . . . , rrtx) of vertices. We can think of this as a "time evolution" of the quiver, and 
for our case at hand will be related to the geometry of Figure 2. We define the quiver at 
"time" t by 



Q(t) := p.mtP-mt-1 ■ ■ ■ Arn-i<3 , Q{0) := Q . (2.10) 

cluster 
(Q,m) 



We can then define the cluster partition function z 7 ^? by 



^«S:=HA mi ...A m Jout>, (2.11) 

for the initial and final states |in) £ "Hq(o) and |out) £ Hq(l)- The partition function depends 
on the choice of initial and final states, whose dependency is suppressed from the notation. 

The cluster partition function has been studied in the context of the wall crossing phe- 
nomena of 4d M = 2 theories. The quiver in that context is a BPS quiver 4d M = 2 theories, 
and the partition function (2.11) is the expectation values of the Kontsevich-Soibelman mon- 
odromy operators [40] (e.g. [22]), see also [3, 19]. 



14 In general certain linear combinations of Yj are in the center of the algebra. In the case of the quantum 
Teichmiiller theory discussed in section 2.2, the corresponding parameters specify the holonomies of the fiat 
SL{2) connection at the punctures of the Riemann surface. 
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2.2 Quantum Teichmiiller Theory 

The quantum Teichmiiller theory fits neatly into the general framework of the previous 
subsection [41, 42]. 15 

Let us begin with an ideal triangulation of E, i.e. a triangulation where all the vertices 
are located at the punctures of E. 

Let us choose an ideal triangulation T, i.e., a triangulation such that all the vertices are 
at the punctures. Given a triangulation of a 2d surface, we can associate a quiver by drawing 
a 3-node quiver for each triangle (Figure 5). The index set of this quiver i, j, . . . , 6 / is the 
set of edges, and that the matrix Q satisfies 

Qij €{-2,-1, 0,1, 2}, (2.12) 

for all i,j G I. See Figure 6 for an example. 




Figure 5: Given a triangulation of a 2d surface, we can associate a quiver by drawing the 
quiver in this figure for each triangle. The vertex of the quiver is associated to an edge of 
the triangulation. 
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Figure 6: An example of a quiver defined for a triangulation, for a 5-punctured sphere. 

Given this quiver we could construct the algebra Aq and the Hilbert space T-Lq as in the 
previous subsection. In the Teichmiiller language, y, is the quantization of the Fock (shear) 
coordinate [44], which is the coordinate of the Teichmiiller space. The commutation relation 

15 The connection of TBA and 3-manifolds has been pointed out in [43], long before cluster algebra was 
formulated. 
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(2.5) represents the standard symplectic form (Weil-Petersson form) on the Teichmiiller 
space, and the Hilbert space %q coincides with the Hilbert space of the quantum Teichmiiller 
theory. 

The description to this point relies on the choice of a triangulation. We can change the 
triangulation by a flip, a change of a diagonal of a square (Figure 7). In fact, it is known 
that any two ideal triangulations are related by a sequence of flips. 



4 4' 




Figure 7: (a) a flip changes a diagonal of a square, consisting of two triangles in the trian- 
gulation. (b) a flip corresponds to a mutation of the associated quiver. 



It is easy to see that the effect of such a flip is translated into a mutation of the associated 
quiver diagram. In this case, the mutation rule (2.8) simplifies to 

yi yi(l + q y 5 ) , y 2 -»■ y 2 (l + (f 1 ^ , ^ _ 

ys -> ys(i + g ys) , y4 -> y^i + g^Vs > ys -> y^ 1 , 

where we use the labeling in Figure 7. 

As discussed in introduction, we consider the effect of an element of the mapping class 
group ip on "Hg. This changes the triangulation, which in turn is represented by a sequence 
of flips m = (mi, . . . , m L ). We could then define the associated partition function Z?-^? 1 to 
be Zf$% defined in (2.11). 

Note that given </? the choice of flips m is far from unique. However, different choices of 
m for a given ip leads to the same partition function, thanks to the quantum dilogarithm 
identities [39, 45]. 16 



16 



Wc need to keep track of the labeling of edges in order to write down quantum dilogarithm identifies. 
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2.3 Cluster Partition Functions 



Let us next evaluate the partition function (2.11). In order to convert the operator product 
in (2.11) into numbers, we need to insert complete set of basis in between the operators, as 
in standard in quantum mechanics. Namely we choose a polarization in Hq, i.e. choose a 
set of coordinates x and momenta p. Then the Hq has a standard representation on the 
coordinate/momentum basis \x), \p), and by inserting a complete basis 

dx \x){x\ , (2.14) 

and thereby converting operators into c-numbers. 

This was precisely the analysis of [13], which discussed the once-punctured torus bundles. 
However, in that paper we need to insert difference complete sets with different polarizations 
((2.14) with different x's) depending on the order of operator products, and a careful case- 
by-case of analysis was necessary to take care of this. In this section we generalize this 
analysis to general mapping cylinder in a more uniform manner, drawing heavily on the 
results from [45] in the context of quantum/classical dilogarithm identities. See also [14, 46] 
for the semiclassical counterpart of our analysis. 

To choose a canonical choice of polarization in Hq, let us prepare a set of variables u,, p; 
for all the edges i, satisfying commutations relations [u$, Pj] = \f—\^b 2 b it y Define y,; by 



y< : = ex P ^Pi + ^2Qij u ^j 



(2.15) 



It then follows that y i 's satisfy the commutation relation in (2.5), justifying the notation. 
The variables Uj, Pi have standard representations on the basis \u), \p): 

\Xi\u) = Ui\u) , pi\u) = -V^l7rb 2 — \u) , 

,- a ( i \ < 2 ' 16 > 

Ui\p) = V^lnb 2 — \p) , Pi\p) = Pi\p) , (u\p) = exp up , 

dpi \y/-lnb 2 J 

and moreover we have the completeness 

W)H=1, /*W<Pl=l, (2.17) 

with du := dui, dp := Y\ i dpi. 

Note that this is a highly redundant description of the commutation relation (2.5); we 
have doubled the number of variables. However the advantage is that we can canonically 
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evaluate the expectation values of yf 



(u\yi\p) = exp ^pi + ^ Qi,i u jj ( u \p) ■ (2-18) 

After these preparations, we could now evaluate the partition function; the operator fi^ 
now acts in a concrete manner in the states \u), \p), and we could evaluate its expectation 
value by inserting the complete sets. We can go through this this exercise following [45], 
see also [3]. 17 18 The answer depends on the choice of initial and final states. We here take 
these to be in the u-basis: 

| in) = |u(0)) , | out) = \u{L)) . (2.19) 

We can evaluate the partition function in different initial and final states, by applying con- 
verting the expression to the basis of (2.19): 

(in|/} mi . . . [i mh | out) = J du(0)du(L) (in|u(0)) («(0) \fi mi . . . fi mL \u(L)) (u(L) |out> . (2.20) 

We will come back to the change of initial and final states in section 4.2. In the basis of 
(2.19), the answer reads 



L-1 L-1 L-1 



z IqS= I ri n II * ( irb\ PmM ~^ QnH ' j{t)Uj{t) 

J t=0 t=l t=o \ * 



(2.21) 



( -^(l 



where 



ft(f + l):=< - PmM il = mh (2.22) 

+ [QmtA^+Pmtit) (l ^ m t) , 

and eft(^) is the quantum dilogarithm function defined in Appendix A. For notational sim- 
plicity we did not explicitly show some of the indices I; for example u(t)p(t) : = 

The expression (2.21) has a large number of integral variables Ui(t),pi(t). However, most 



17 In this evaluation we decompose the action of fik into two parts, a linear exchange of variables and 
a conjugation by a quantum dilogarithm, see [39, 45, 13]. In the language of quantum mechanics, this is 
the translation from the Heiscnberg picture to the Schrodingcr picture. The quantum dilogarithm in (2.21) 
comes from the operator representing this conjugation. 

18 We need to modify the argument of [45] slightly to incorporate in and out states. We do not need to 
include their v* , which represents the re- labeling of the edges and is crucial for the quantum dilogarithm 
identifies of [45] but not for the purposes of this paper. Note also that $>b(%) in [45] is our eb(x)^ 1 . 
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of them can be trivially integrated out. Indeed, the power of the exponent in (2.21) reads 
u(t)p(t) - u(t + l)p(t + 1) = u mt (t)p mt (t) + u i(t)Pi{t) 

+ U mt+ i(t + l)p mt (t) - ^ + + IQrntAW+PmM) ■ 

In particular the variable Pi(t) (i ^ m t ) does not appear inside the argument of dilogarithm, 
and only appears in the linear term (v,i(t) — Ui(t + l))pi(t). Integrating out pi(t) (i ^ m t ), 
we have 



Uiit) = Ui(t + 1) (i ^ m t ) 
and hence most of the u- variables are identified, leading to 



(2.23) 



^cluster 
Z (Q,m) 



L-l 



L-l L-l 



-I 



ldPmMl[du(t)\ 



e 6 



t=o 



X 



t=l 



t=0 



2nb 



Pn 



exp ( ^^p mt (t)(u mt {t) +U mt (t + 1) - ^[Qm t ,i(*)]+ U <(*)) 



(2.24) 



We can integrate over p mt (t) by (A. 7) in Appendix, leading to (up to a constant overall 
phases irrelevant for the identification of 3d M = 2 theories) 



^cluster 



iidu(t) n^(^^)-f) 'ex P 

t=i t=o v ' 



(2tt&) : 



Z'(t)Z"(t) , (2.25) 



where we defined Z'(t), Z"(t) by 



Z\t) := 2 
:= 2 



i 

(t) + M mt (t + 1) - Y^[-Qm t ,i( t )]+ U i( t ) 



(2.26) 



and 



(2.27) 
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For later purposes we also define 



Z(t) := V^lKbQ - 2 



= V^lTrbQ - 2 Qm*,i(*H(*) J , (2.28) 

satisfying 

Z(t) + Z'(t) + Z"{t) = yf—litbQ, . (2.29) 



The factor ^J^-i^lbQ will turn out to be useful for later considerations. 

In (2.25), we need to remember that the variables Ui(t) are constrained by (2.23). To take 
these constraints into to account, it is useful to introduce a notion of a mutation network. 

2.4 Mutation Networks 

Let us again begin with a pair (Q,m), a quiver Q, and sequence of mutations m = 
{m 1 ,m 2 , ...,m L }. We then associate a graph, a mutation network (see Figure 8). 

The graph is always bipartite, i.e., vertices are colored either black or white and edges 
connect vertices of different colors. We denote the set of black (white) vertices of the network 
by B (W). 

A black vertex represents a mutation, one of the m's. A white vertex represents integral 
variables itj(t) of the previous subsection. 

The black vertex m G B, representing a mutation m, is connected to two white vertices 
(denoted by x^ in \x'^ rn ' in Figure 8) by dotted lines — these represent the variables u mt (t) 
and u mt (t + 1) in the previous subsection. A black vertex is also connected with other white 
vertices by undotted lines — these white vertices represent Uj(£)'s with i ^ m t . The number 
of undotted arrows from m G B to w 6 W is determined by Q m ,w '■= Q(t) muW , which is one 
of the components of the quiver adjacency matrix before the mutation. In general there are 
multiple arrows from m to w (or from w to m, depending on the sign of Q m ,w)- 

This rule defines the mutation network. Around a black vertex, the network represents 
the mutation and the quiver vertices affected by it. Around a white vertex, the network 
describes the creation of a new integral at some time, and its annihilation at a later time. 
The network is naturally concatenated when we combine two mutation sequences. Concrete 
examples of mutation networks will appear later in section 3.4. 

Let us come back to our partition function (2.25). We learn from (2.23) that the inde- 
pendent variables are associated with the white vertices of the mutation network, we denote 
this variable by x w (w G W). Some of the edges are in the in the initial (final) quiver Q(0) 
(Q(L)), and others not; we denote the difference by w G W exk if they are, and w G W int if 
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b 




x >(m) 

Figure 8: (left) A mutation network around a mutation m. The mutation is represented by 
a black vertex m, and the integral variables are represented by white vertices. Two of the 
white vertices, x^ m ' and x ,( - m \ represent the variables associated with the mutated vertices, 
and are connected by dotted lines. Other white vertices are connected by undotted lines, and 
comes with the charges Q m , w - (right) When the quiver arises from an ideal triangulation of 
a 2d surface (section 2.2), a mutation network always looks as in this figure around a black 
vertex. 

they are not. By definition we have 



w cxt u w int = w, w ext n w int 

We can now rewrite the result (2.25) as 



^cluster 
<(Q,m) 



dx w [ e b f 7r - r Z'{m 



-1Q 



2nb 



cxp 



(2.30) 



where we defined 



Z'(m) := 2 
Z"(m) := 2 



. x ( m ) _ x '( m ) 



^ ^ [Qm,'i 



x {m) + x /(m) 



£[-« 



(2.31; 



and 



Z(m) := V^lirbQ - 2 



Eft' 



m,w\ +%w 
,w£W w£W 



-1tt6Q-2 I Qm, w x w \(2 
\wew ) 



Note that the integrand of (2.30) factorizes into contributions from each mutation (m). The 
expression (2.30) will be the crucial ingredient for the identification of 3d gauge theories in 
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section 4.2. 



When we discuss 3-manifolds (and associated 3d M = 2 theories), we concentrate on the 
case when the quiver is determined from an ideal triangulation of a Riemann surface (section 
2.2). In this case the mutation network always looks as in the right of Figure 8 around a 
white vertex, namely two lines corresponding to mutated vertices, and four lines with two 
charge +1 (denoted by a^ m \ c^ m ') and two —1 (denoted by b^ m \S m ^). 19 In this notation we 
have 

Z(m) = V^lnbQ + 2 (-a (m) + c (m) + 6 (m) + S m) ) , Z'(m) = 2 (-x (m) - x' (m) + a (m) + c (m) ) , 
Z"(m) = 2 (x (m) + x /(m) - b {m) - d {m) ) , (2.33) 

where for notational simplicity we used a*-" 1 -* also for the associated variable, which should 
be written x a ( m ) in our previous notation. Interestingly, this (in particular the variable Z(m) 
is precisely the coordinate transformation from cluster x-variables to the cluster y-variables, 
or in the Teichmuller language from the Penner coordinates (geodesic length) to the Fock 
coordinates (shear coordinates, i.e. cross ratios), see Table l. 20 




Jirn) c (m) S m) c (m) 



Figure 9: Mutation of a quiver at vertex x, representing a flip of an 2d ideal triangulation. 



Table 1: Dictionary between mutation networks, Teichmuller theory and cluster algebras. 
The Teichmuller interpretation does not exist for general quivers which does not arise from 
an ideal triangulation of a 2d surface. 



mutation network 


cluster algebra 


Teichmuller space 


variable x w (w G W) 


cluster x-variable 


Penner coordinate 


variable Z(m) (m e B) 


cluster y-variable 


Fock coordinate 



Coming back to the general expression (2.30), we find that the integrand depend on 

19 a( m ) and cS m \ or 6' m - ) and S m \ could be identified with each other, as in the case of the quiver coming 
from the triangulation of the once-punctured torus. 

20 A numerical factor of 2 in (2.33) is consistent with [45, Remark 5.1]. 
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2\B\ variables Z'(m), Z"(m). As we mutate the quiver, the number of such variables grows 
roughly twice the speed as the number of variables x w 's (|Wi nt |), since we obtain two new 
variables Z'(m),Z"(m) for each mutation. 21 This means that there should be constraints 
along Z'(m), Z"(m)'s. 

In fact, we find the following constraint for each w £ Wi nt , which eliminates the number 
of variables by | W int | ! 

Suppose that we fix a white internal vertex w £ W- m f This vertex is connected with 
many mutations m £ B, and two of mutations are connected with w by dotted lines — a 
vertex is created by a mutation, and deleted by a mutation after some steps. We call these 
mutations to be of type 1. The remaining mutations are called of type 2 (of type type 3) if 
Q m ,w < (if Qm,w >0). We then find 

Z{m)+ [-Qm, w }+Z'(m)+ [Qm,w]+Z"(m) = 2nbQV^l . (2.35) 

■m: typcl m: type2 m: type3 

This relation will correspond to the superpotential constraints for R-charges in 3d M = 2 
theories in section 4.2. 

To show (2.35), it is useful to go back to the notation of (2.25), under which (2.35) reads 

Z(tx) + Z(t 2 ) + Yl 1-Qmj + Z'(t) + [QmJ + Z"(t) = 2nbQV^l , (2.36) 

ti<t<t 2 h<t<t 2 

where we used the notation that the variable Ui(t), corresponding to w, is generated at t = t\ 
and annihilated at t — t%, i.e., i = m tl = m t2 . Using the definitions (2.26), (2.28), the left 
hand side of (2.36) becomes 

mt,i \""mt 

(t) + U mt (t + 1)) 

v 3 j J h<t<t 2 

+ ([QithMQitUjU - [Q(t)3,m t ]+[Q(t)mJ+) + 7r&gV=l 

tl<t<t 2 j 

(2.37) 

From the definition of mutation (2.2) and the constraints on the variables Ui(t) (2.23), the 
third term inside the bracket of (2.37) is equivalent with 

Y + ^ - u ^ = E E + ViMtt + !) - • 

ti<t<t2 j y^rnt ti<t<t2 jj^mt 



21 In general we have the relation 

|W int | + |W ext | = |W| = |7| + |S| . (2.34) 
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The sum over j becomes over all j & I when this is combined with the second term in (2.37). 
After many cancellations, this leads to 

(Q&kjUjih) - Q(h + + 1)) = (QfahMfo) + QitihjUjih)) . 

j j 

This cancels the first term in (2.37), with the only remaining term in (2.37) being the fourth 
term, giving 2irbQ\/—l. This proves (2.35). 

3 3-manifolds 

In this section we discuss in more detail the geometry of hyperbolic 3-manifolds determined 
from the action of the mapping class group on quantum Teichmiiller space. We have already 
seen in section 1 that the the geometry is a mapping cylinder or torus, however we wish to 
obtain more detailed information of the geometry, to make contact with the discussion of 
section 2. 

The goal of this section is threefold. We first discuss the canonical ideal triangulation 
of our 3-manifold (section 3.1), which originates from an ideal triangulation of a 2d surface. 
Second, we discuss how to modify the construction to obtain more general geometries, by 
identifying unglued faces of a mapping cylinder (section 3.2). Third, we discuss the gluing 
equations of hyperbolic tetrahedra, and show that it arises from the semiclassical limit of 
the partition function discussed in 2 (section 3.3). 

3.1 Canonical Ideal Triangulations 

Let us discuss the ideal triangulation of M. For a mapping cylinder there is a canonical 
choice of ideal triangulation [47, 48, 49]. 

As we have seen already, the action of tp could be traded for a sequence of flips (which 
in turn is identified with a mutation sequence m). We can then associate a tetrahedron 
for each flip; given a 2-manifold with triangulation, we can attach a tetrahedron (squeezed 
like a pillowcase) and we effectively obtain a new 2-manifold with a different triangulation, 
related with the original one by a flip (Figure 10). By repeating this procedure we obtain a 
sequence of tetrahedra, whose faces are glued together. The 3-manifold is now decomposed 
into tetrahedra: 

M = |jA m . (3.1) 

m 

Our mutation network contains all the information about canonical triangulations. We 
associate an ideal tetrahedron for each mutation m inside m (represented by a black vertex). 
Since a tetrahedron has six edges, each black vertex is connected with six white vertices (see 
the right figure of Figure 8). The mutation network also specifies how to glue these tetrahedra 
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together, and hence the gluing equations in Appendix B; an edge w (a white vertex) is shared 
by all the tetrahedra (black vertices) which are connected with the w in the mutation network 
(see also Figure 18 in section 3.3). 




Figure 10: Attaching a 3d tetrahedron to a 2d triangulation flips the diagonal of a square. 

The canonical triangulation is so far a combinatorial triangulation, however, we can to 
promote it to an ideal triangulation by the hyperbolic tetrahedron when the 3-manifold is 
hyperbolic. 22 This means that each tetrahedron is an ideal tetrahedron in H 3 , and there is 
a (complete) hyperbolic structure of the 3-manifold (see Appendix B for brief summary of 
3d hyperbolic geometry needed for this paper). 

The mutation network in the 3-manifold cases discussed in this section is reminiscent of 
the "braid/tangle" of [3]. This is the branched locus of the IR geometry, which is a double 
cover of our 3-manifold M. In both cases we associate a basic building block, either a black 
vertex (for mutation network) or a crossing, to a tetrahedron (Figure 11). 

However, it is important to keep in mind that the "braid" in their paper, or rather 
the branched locus, is not the braid/knot discussed in our paper. In fact, inside an ideal 
hyperbolic tetrahedron our knot (which appear in "knot complement"), for example, goes 
through the vertices of tetrahedra, whereas the "braid" in [3] goes though the faces of 
tetrahedra (Figure 11). In other words their "braid" path through the zeros of quadratic 
differential of a 2d surface in the section of the 3-manifold, whereas our braids pass through 
the poles. In the rest of this paper the words "knot/link/braid" will alway refer to the 
knot/link/braid in our sense. 

3.2 Capping the Braids 

There is a caveat in the discussion to this point. The 3-manifold obtained in this construction 
is of special type S x J, and does not seem to be general enough. However, what saves the 
day is that by suitably identifying boundaries of this 3-manifold it is possible to obtain a 
rather large class of 3-manifolds, including all the link complements in S 3 . 

What we do here is to identify unglued faces of tetrahedra on the boundary of the mapping 
cylinder. Depending on the identification we obtain different 3-manifolds. 

22 A 3-manifold is "generically" hyperbolic; a knot complement in S 3 , to be discussed in next subsection, 
is hyperbolic unless the knot is one of the torus knots or their satellites. 
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flip tetrahedron branched mutation network 




Figure 11: A flip corresponds to a tetrahedron, which is represented as a crossing in the 
branched locus of [3], or a black vertex of the mutation network introduced in this paper. 
Beware of the difference between the branched locus of [3] and the braids in this paper. 

Such an identification has been worked out for the case of the 4-punctured sphere [50, 51]. 
In this case, there are four faces in the triangulation, and we first identify two of the faces 
and then the remaining two (see Figure 12). We can verify that this face identification gives 
rise to the identification of the braids passing through the four punctures, and hence the 
braids are capped off into links. We can realize a large class of knots called 2-bridge knot in 
this way. 

We can generalize this construction to 2n-punctured spheres (Figure 13). In this case 
we can again identify the faces of the boundary surface, leading the the identification of the 
braids. By applying an element of the mapping class group we could obtain an arbitrary 
identification of 2n-braids. 

In particular we could choose the identification as in Figure 14 for the 2n-punctured 
sphere bundle, both for the in and out states. A link obtained after such an identification 
is said to have 2n-plat representation. 23 It is known that an arbitrary link has a 2n-plat 
representation for some n [52], which means that our procedure include arbitrary link in S* 3 . 
It is clear from the Figure 14 that the resulting link is determined from an element of the 
braid group £>2„ (recall also Figure 3). 

The identification of the faces induces identification of edges, which is translated into 
the identification of corresponding variables in Aq and Hq. 24 In mutation networks 
this procedure of gluing boundary faces of mapping cylinders is simply translated into the 
identifications of white vertices (edges of tetrahedra). We will discuss the example of the 
figure eight knot complement in section 3.4. 

23 This is similar, but different from the so-called braid representation of a knot/link. 

24 Strictly speaking this identification could involve factors of q as long as they become trivial in the 
classical limit. These factors affect our partition function, but not the semiclassical analysis of section 3.3. 
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Figure 12: By identifying the faces of 4-punctured sphere bundles, we can identify the four 
braids (corresponding to four punctures) in pairs; the same figures can be found in [50, 51]. 




Figure 13: Similar identification exists for 2n-punctured sphere bundles, by iteratively mak- 
ing a pair among the faces. This is an example of n = 3. 
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Figure 14: A 2n-plat representation of a knot/link. The link is determined by an element 
of the the braid group B 2n - The example on the right is determined by an element sfs^sf of 
Bq, and is a link often called 61. 

Comparison with Heegaard Decomposition 25 

Our capping procedure is closely related with the Heegaard decomposition of a 3-manifold, 
and its generalization. 

A Heegaard decomposition states that a closed 3-manifold M has a decomposition of the 
form 

M = H 1 U„H 2 , dH x = dH 2 = S,, , (3.2) 
where H\ and H 2 are the handlebodies 26 and ip is an element of the mapping class group of 

The decomposition (3.2) could be represented as in Figure 16 (a). As the figure shows, 
the only effect of the handlebody should be to choose a specific element 

|handlebody) G Hz , 

and we could evaluate the partition function by substituting these states in the in and out 
states (1.5). 

For our purposes we yet need a small modification; we need to include knots, and consider 
a 3-manifold with torus boundaries. In the two-dimensional slice the knots are the point-like 
in two-dimensional surface S, and serve as a puncture of £ (recall Figure 3). 

25 This part is outside the main track of this section and could be skipped on first reading. 
26 Colloquially they are the "simplest" 3-manifolds with boundary S (Figure 15 (a)). For example, the 
handlebody for the two-sphere S 2 is the three-dimensional ball B 3 . 
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Figure 15: A handlebody (left) and a tanglebody (right). 

Correspondingly we need to consider handlebody, with knots inside. We call these a 
tanglebody, see Figure 15 (b). 27 Note that tanglebody exists only when the number of 
punctures of £ is even, since whenever a knot comes into the tanglebody it need to come 
out. It is also clear that given E the tanglebody is not unique. For example, in the case of 
the 4-punctured sphere we have three tanglebodies of Figure 17, and each of them give rise 
to different states. Note that corresponding choice is present in Figure 12, when we identify 
the four faces of 4-punctured spheres in pairs. 




M 

Figure 16: Heegaard decomposition (a) and tanglebody decomposition (b). 

27 As commented already, a handlebody for a sphere is simply a ball, and in this case it is straightforward 
to define the corresponding tanglebody as a ball with braids deleted from there. 
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Figure 17: Three tanglebodies for the 4-punctured sphere. 



Now we can generalize the Heegaard decomposition and consider the "tanglebody de- 
composition" (Figure 16) 



where T 12 are tanglebodies. 

Our gluing procedure explained above (Figures 12, 13) should be essentially the same 
as gluing the tanglebodies, in the sense that in both cases the braids are identified in the 
same way. The generality of the Heegaard decomposition roughly explains the generality of 
our approach. It would be interesting, however, to understand the relation between the two 
approaches in more detail. 

3.3 The Semiclassical Limit 

Finally, let us directly verify that the partition function discussed in the previous section 
reproduces the gluing equations (B.4) of the associated hyperbolic 3-manifold. This is the 
direct demonstration of the consistency between this section and section 2. The semiclassical 
analysis of our partition functions can be found in [45], 28 and the fact that the saddle point 
analysis of our partition function reproduces the gluing equations, as well as the connection 
with cluster algebras, has already been worked out in [14], see also [46]. 

The classical limit of the Chern-Simons theory is t — ?• oo, or equivalently b — > 0, q — > 1 
limit of (1.1). It is straightforward to take the b — > limit of the expression (2.30) obtained 
in the previous section, and we find (using (A. 5)) 



where W(x w ) is written as a sum over the contributions W^(x w ), each associated with a 

28 Our semiclassical analysis here is actually slightly different from [45] in that we have used (2.25) with 
the variables pi(t) integrated out, whereas they used (2.21) and extremized also with respect to Pi(t). Both 
methods give essentially the same results. In fact, the proof of (2.35) here is somewhat similar the proof in 
[14], although proven in different variables. 



M = 7\ U v T 2 , 



(3.3) 




(3.4) 
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flip m: 

= E ( Li 2(^' M ) " \z\m)Z\m) \ . (3.5) 
The saddle point of this integral is given by 

expf— J=l, weW int . (3.6) 

We now claim that this equation is identical to the gluing equation for the canonical ideal 
triangulation, and that W(x w ) is a generating function of the gluing equations described in 
[33]. 

Let us pick up a particular edge of the triangulation. This is represented by a vertex 
vertex w G B. Then the x^-dependent part of W is a sum over contributions (denoted by 
yy( m )) from all the tetrahedra m containing the edge e, where in this case is given by 

(3.5). 

The mutation m is divided into three types, of type 1, type 2 or type 3, as discussed in 
section 2.4. The x^-derivative of the W*-"^, in each of the three cases, are given by 

21og(^ m >) - i7T - Z(m), [-Q m , a ]+(21og(^ m >) - Z'(m)), [Q m ,J + (2 log(^) - Z"(m)) (3.7) 

for type 1, type 2, and type 3, respectively, where we introduced tetrahedron modulus for 
the tetrahedron m by 

and we introduced the three parameters z,z', z" related by 

z >(m) = 1/(1 - ^)), Z " {m) = 1 - Z^- 1 . 

These will correspond to three different parametrization of a tetrahedron, as explained in 
Appendix B. 

When we collecting these factors and sum over m, the terms linear in Z(m), Z'(m), Z"{m) 
in (3.7) cancel out due to (2.35) (note bQ — > 1 in the semiclassical limit). This means that 
we are left with 

exp (^^\ = n z ^ n (z'M)hQ m ,j + n ^^q^ . (3 . 9 ) 

^ / m: type 1 m: type 2 m: type 3 

When our theory is associated with the 3-manifolds, this is exactly the gluing equation 
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(B.4) in section 3.1, and our derivation represents the fact (proven in [14]) that y- variables 
automatically solve the gluing equations. 

Interestingly, there is a natural symmetry cyclically exchanging the Z(m), Z'im) and 
Z"(m). In fact, one can replace (3.5) by 

WW = ( Li 2(e Z(m) ) - \{Z{m) - V^ln)Z\m)\ , (3.10) 

and we can verify that this still gives the correct saddle point equation. For the 3-manifold 
cases we can understand this symmetry as the cyclic exchange of three modulus parametriza- 
tion z,z',z" ((B.3) in Appendix B). 

The gluing equation has a rather concise expression in the mutation network defined in 
section 2.4. The gluing equation can be written down for an internal edge of the triangulation, 
and hence is associated with a white internal vertex of the network, see the left of Figure 
18. The part of the mutation network around a white vertex is in direct correspondence 
with the projection of shape of tetrahedra along the corresponding edge (Figure 18), or 
equivalently the boundary torus around the edge. Note that mutation network also specifies 
the parametrization of tetrahedra, i.e., whether we use z, z' or z" (Figure 19). 




Figure 18: The mutation network around a white vertex represents how the tetrahedra are 
glued together. 

3.4 Examples 

For concreteness and for comparison with earlier results, we here work out two 3-manifold 
examples. This will illustrate the generality and usefulness of our approach. 

3.4.1 Figure-Eight Knot Complement 

Let us first discuss one of the most famous hyperbolic knot complements, the figure eight 
knot complement. This can actually be realized as a mapping cylinder, and is discussed in 
detail in the paper [13]. We here use the 4-plat representation of the knot. 
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Figure 19: The mutation network specifies the figure of the boundary torus, as well as the 
parametrization of the shape parameters of the ideal tetrahedra. The three different types 
of edges correspond to three different parametrization of a tetrahedron modulus. 



In the standard 4-plat representation of the figure eight knot, we need four generators 
of braid group B4. For practical computations, however, it is more efficient to incorporate 
some of the mapping class group actions into the choice of the caps (recall section 3.2). We 
can then realize our knot by a single flip with caps on both ends, leading to the mutation 
network in Figure 20. Interestingly, this gives rise to the famous ideal triangulation of the 
figure eight knot complement by the two ideal tetrahedra, found in [53] (Figure 21): 

A = A', B — B', C = C, D = D' . 

with the identification of edges: 

x :— a — d — f — a' — d' — f , y := b = c = e = b' = c' = e' . 

The semiclassical limit of the partition function (3.5) gives 

W(x, y) = Li(e 2( - e - e ' +fc+d) ) - ^ 2(-e - e' + b + d) 2(e + e' - c - /) 

+ u(e 2{ - a - a ' +h+d) ) - i 2(-a -a' + b + d)2(a + a' - c - /) 
= 2(Li(e 2x - 2y ) + 2(x - yf) , (3.11) 

where x,y are parameters associated with the two edges. When we define z := e 2x ~ 2y , the 
critical points are given by 



z — z + 1 = , i.e., z 



-1 ± 



This corresponds the complete hyperbolic structure (and its complex conjugate) of the figure- 
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Figure 20: The mutation sequence, the mutation network and the knot projection for the 
figure eight knot complement 4^. This gives the ideal triangulation in Figure 21. 
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gluing 




eight knot complement. We can verify that the critical value of W reproduce the hyperbolic 
volume and the Chern-Simons invariant of the 3-manifold. 

By following similar methods, we can compute the partition function for any link com- 
plements in S 3 . See Figure 22 for another example. The general recipe for reading off a 
mutation sequence for a given Dehn twist is explained in Appendix C. This rule, in our 
experience, is rather useful for practical computations. 



3.4.2 Once-Punctured Torus Bundles Revisited 

To illustrate the usefulness of our formalism, let us work out the example of the once- 
punctured torus bundle. This example has been worked out in detail in [13]. In particular it 
was found there that the quadratic piece of the W depends in a subtle way on the mutation 
sequence. We here find that the rules proposed in this paper reproduce the findings in [13]. 

Let us discuss the mapping cylinder for the once-punctured torus with ip = LLRRRRL, 
in the notation of [13]; this is sufficient to discuss the general pattern. The mutation network 
is given in Figure 23. We use the form of the semiclassical potential in (3.10) 



t=o 



Li 2 (e 



Z{t)y 



;(Z(t) 



br)Z'(t) 



(3.12) 



with 



Z(0) - A 


/-Ttt 


= 2b- 


-2c, 


Z'(0) = - 


-a — a' + 2c, 


Z(1)-a 


/-LtT 


= 2b- 


- 2a, 


Z'(l) = - 


-c - d + 2a', 


Z(2) - a 


f-\-K 


= 2d 


-2a', 


Z\2) = 


-b-b' + 2a, 


Z(3) - a 


f-i-n 


= 2b' 


-2a, 


Z'(3) = 


— c — c + 2a , 


Z(4) - a 


/-Ltt 


= 2c" 


- 2a, 


Z'(4) = 


-6" - 6"' + 2a', 


Z(5) - a 




= 2b" 


- 2a', 


Z'(5) = 


-c" - c'" + 2a', 


Z(6) - a 


/-TtT 


= 2b" 


- 2c"' 


Z'(6) = 


: -a" - a' + 2c'" 



(3.13) 
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Figure 22: The mutation sequence, the mutation network, knot projection and the boundary 
torus for the knot complement 52- Compare this with [51, Appendix]. 
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where we have again used the simplified notation w for the variable x w , and for mapping 
torus we identify a" = a, b" = 5, d" = c. As have already seen, the saddle point of this 
potential W gives the gluing equation of the hyperbolic 3-manifold. 



a 

O 




Figure 23: The quiver, the mutation network and the boundary torus for a once-punctured 
torus bundle with ip = LLRRRRL. Compare the boundary torus with those in [49, 13]. 

Let us define w t := Z(t — 1) — y/^lir for t = 1, . . . , 7, and re-express W in term of Wt's. 
We then find 

W{w) =Li 2 {-e Wl ) + -w x (-w 7 -w 2 + 2 Wl ) 



+ Li 2 ( 


-e W2 ) 


+ ■^w 2 (-w 1 + w 2 + w 3 ) 


+ Li 2 ( 


-e W3 ) 


+ ^^3(^2 + ^4) 


+ Li 2 ( 


-e WA ) 


1 / \ 

+ -W4,{W 3 + W 5 ) 


+ Li 2 ( 


-e W5 ) 


1 / 

+ ^^5(^4 + we) 


+ Li 2 ( 


-e m ) 


1 / \ 

+ -W 6 (^5 ~W 7 + Wq) 


+ Li 2 ( 


-e W7 ) 


1 . 

+ -w 7 (-w 6 - w 1 + 2w 7 



Note that associated with the Li 2 (— e Wt ), there is a quadratic expression with respect to 
Wt-i, Wt, Wt+i- This is determined by the choice of either L or R for the neighboring (t — 
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1), t, (t + l)-th flips, which for t — 1, . . . , 7 are given by 



1 : LLL , 2 : LLR , 3 : LRR , 4 : RRR , 5 : RRR , 6 : RRL , 7 : RLL . 

We can compare this with the results of [13], and find complete agreement. For example, for 
the first line of (3.14) we have LLL, and this coincides with the case 1 of [13, section 3.4]. 

As this examples illustrates, while it is possible to write the final expression only in 
terms of cluster y-variable, the most systematic expression requires the use of the cluster 
x-variables. 



In this section we construct 3d = 2 theories T[(Q, m)] satisfying (1.2), based on the 
results of section 2. 

4.1 Sf Partition Functions 

Let us first summarize the basic ingredients of 3d M = 2 theories, and their partition 
functions [9, 10, 11, 12]. 

Our theories have a number of U(l) symmetries, which are labeled by the set H. Some 
of the U(l) symmetries are flavor symmetries, and others gauge symmetries. We denote this 
by i G F and i G G, respectively. By definition we have FUG = H , and Fn G = 0. For each 
i there is a corresponding parameter a iy the scalar of the associated Af = 2 vectormultiplet. 
When i 6 F, the correspondingly vector multiplet is non-dynamical and the parameter <7j is 
called a real mass parameter. 

We also include Chern-Simons terms, including off-diagonal ones. This could be described 
by a symmetric matrix kij for i,j G H: 



for dynamical/background gauge field A^j. The Chern-Simons term has to obey a quanti- 
zation condition for the invariance under the large coordinate transformation and for the 
absence of parity anomaly, and in particular fcjj should be a half-integer. Note also that 
we consider fcjj for either % or j in F; these are Chern-Simons terms for background Chern- 
Simons terms, which plays crucial roles when we gauge the associated global symmetry. 

The partition function of a 3d M = 2 theory depends on the Chern-Simons term 
ki j G H) and the real mass parameters cr^ (i G F). Here is a 1-parameter deformation 
of the S* 3 preserving U(l) x U(l) isometry. More explicitly it is defined by [12] 



4 3d AT = 2 Theories 




(4.1) 



b 2 \ Zl \ 2 + b- 2 \z 2 \ = r 2 , 



(4.2) 
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with Z\,Z2 G C. The partition function turns out to be independent of r. 

Now the S* 3 partition function, after localization computation, takes the following form: 

% [Kj (i,jeH),&i (i£F)] = J ^Y\.d(Ti^ ^classical ( k, a) Zi_i oop (cr) . (4.3) 

The rules are summarized as follows (we specialize to Abelian gauge theories in this paper): 

• The integral is over all the Abelian global symmetries o"j, (i G G). 

• The classical contribution Z c i assica i(/c, a) is determined by the Chern-Simons term, 29 
and is given by 30 



^classical (k,a) = exp I -y/^ln ^ I 



(4.4) 



The 1-loop determinant Zi_i oop has contributions only from M = 2 hypermultiplets. 
This is given by 

Zl-loop(o-) = Y\ S b (y^,Qi, m O~i + - Qm) ) , (4-5) 

ui: hyper \i(HH / 

where we assumed the Af = 2 hypermultiplet m has charges Qi, m under the U(l)i 
symmetry, and has R-charge q m . The correct value of the IR superconformal R-charge 
is dependent on the mixing of the UV U(l) R-symmetry with flavor symmetries 31 , 
and we can write (4.5) as the holomorphic combination 32 



2 

m: hyper \ iS-ff 



with 



x Q 

Re(cTi) = <Ti, Q i)Tn Im(cTj) = — q m . (4.7) 

Note that in the partition function (4.3) the only distinction between <j{ for % G F and 
those for z G G is that there is we do not integrate over the former. This means that effect 



29 There are also classical contributions from FI parameters. 
kij here is the bare Chern-Simons term, not the effective Chern-Simons term obtained after integrating 
out massive matters. 

31 The correct mixing is determined by -F-maximization [10]. 

32 The reason why the Coulomb branch parameter (real mass parameter) and the R-charge appear in a 
holomorphic combination has been clarified in [54] in the background supcrgravity formalism. 
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of gauging a global symmetry is simply to integrate over the corresponding <7; in the S% 
partition function. 



4.2 Properties of 3d N = 2 Theories 

Let us finally comment on the properties of our theories T[(Q, m)]. We will content ourselves 
with the comments on the basic properties and defer the detailed analysis of our theories to 
a future publication. 

The theory T[(Q,m)} in this paper is defined by the relation (1.2), by comparing the 
partition function (4.3) with our partition function (2.30), which we rewrite here to be 
(using (A. 3) and (A.4)) 

Z {qZ) = I H dw H S b (-^Y^ - °rn) eX P ( f^T^ ~ °™) ~ ^™™°'m 

(4.8) 

where we defined 

and we have again neglected overall phase factors. 

The properties of the theory T[(Q,m)] are summarized as follows. 



1. M = 2 Abelian Vector Multiplets (U(l) Symmetries) 

We associate a U(l) symmetry for each white vertex of the mutation network. This is either 
a global or gauge symmetry depending on whether or not the white vertex is associated with 
the quiver edge in the initial/final states. 

F = W ext , G = W int , H = W. (4.10) 

We also identify the corresponding variables x w and 6^ by the relation 

5, = . (4.11) 

For the case of 3-manifolds, an equivalent way to say this is that the U(l) symmetry is 
associated with an edge of the tetrahedron, and is a global (gauge) symmetry if the edge is 
on the boundary (in the interior) of the 3-manifold. 

Note that not all the U(l) symmetries are really independent — in fact, many of them 
are related by electric-magnetic duality, and this happens when the corresponding variables 
Vj do not commute. Note also that the U(l) symmetries in general have Chern-Simons terms, 
and is determined by the quadratic expression in (4.8). 
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2. J\f = 2 Hypermultiplets 



We associate a M = 2 hypermultiplet $ m for each black vertex m £ 5 of the mutation 
network. The charges Q m>ll) associated with the edges of the mutation network determines 
the charges of this fields under the U(l) symmetries; the parameter a m , which appears inside 
the argument of s&, is a linear combination of <7j's due to the relations (2.26), (4.9) and (4.11). 
In particular the R-charge of $ m is given by 



For the case of 3-manifolds, this hypermultiplet is associated with an ideal tetrahedron. 

3. Superpotential 

Finally, there is a superpotential term among M = 2 hypermultiplet. Given a white vertex 
in the mutation network, we associate a superpotential term. Depending on the 3 types 
(types 1, 2, 3 discussed in section 2.4), the operator O m>w is either a fundamental field or 
involves monopole operators. This means that our theories are generically non-Lagrangian. 
The fact that the superpotential has R-charges 2 could be guaranteed by the relation (2.35) 
we have proven earlier. 

4. The Sp(2N, Z) Action 

Let us come back to the choice of initial and final states of the partition function (2.11); our 
discussion so far assumes the choice (2.19). 

As we discussed in (2.20), the change of the boundary condition induces the change of 
the partition function. For example, when we change (in| from (u(0)\ to (p(0)|, we have 



(p(0)|A mi • • • /4KL)> = / du(0) e ^T u[0)m (u(0)\fi mi . . . f£ L \u(L)) , (4.13) 



which is just a Fourier transformation. More generally, we could can choose a different state 
by an Sp(2N, Z)-transformation, where TV here is given by 2\I\ (Note that we could choose 
to mix variables Ui(0),pi(0) and Ui(L) , Pi(L)) . This Sp(2N,Z) action is lifted to the action 
of the wavefunction, which in turn is identified with Sp(2N, Z) action on 3d M = 2 Abelian 
theories [21]. This involves adding Chern-Simons terms for background gauge fields and 
gauging global symmetries with off-diagonal Chern-Simons terms. 

5. Gauging 

Gluing two 3d M = 2 theories is represented by a concatenation of the mutation network. 
When we have several networks Hi and glue them together at a white vertex w, we gauge 




(4.12) 
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the diagonal subgroup of corresponding global symmetries U(l) 



, and we have 



/ 




(4.14) 



where we showed the dependencies of Zj^ only with respect x w . Such a gauging is necessary, 
for example, when we cap off the braids in section 3.2. 

6. Dimensional Reduction and SUSY Moduli Space 

The semiclassical limit 6—^0 discussed in section 3.3 is the limit where the ellipsoid S* 3 
degenerates into M 2 x Si with a circle of small radius b. Since we take b to be small we are 
effectively reduced to 2d M = (2, 2) theory, but with all the KK modes taken into account. 
The parameters x w will play the role of the vector multiplet scalars (which is complexified 
after dimensional reduction), and W is the effective twisted superpotential obtained by 
integrating out matters. The gluing equation then is the vacuum equation for the 2d theory 
(cf. [55]). The SUSY moduli space of our 3d Af = 2 theory should be constructed from the 
symplectic quotient construction (cf. [56]), imposing (2.35) as constraints. 

7. Loop Operators 

We propose that the variables represent the flavor Wilson/vortex operators for the i-th 
flavor symmetry. In fact, flavor Wilson (and vortex) loop operators wrapping the Hopf fiber 
at the north pole 33 of the base S 2 represent the multiplication and shift to the partition 
function [57, 58], leading to the commutation relation (2.5). This means that insertion of 
the yi in the cluster partition function (2.11) should be identified with the (un-normalized) 
VEV of the corresponding loop operator. The situation is similar to the case of 3d Af = 2 
theories coming from the dimensional reduction of 4d M = 1 theories in [59], which paper 
also proposes to identify (classical) y- variables with the VEVs of loop operators. 

8. Coupling to 4d Af = 2 Theory 

Our quiver in many cases can be identified with a BPS quiver for a 4d Af = 2 theory. 
In these cases it is natural to propose that our 3d theory arises on the boundary of the 
4d Af = 2 theory; the latter couples to the former by gauging global symmetries of the 
former. The BPS wall crossing causes the mutation of the BPS quiver, which is translated 
into the change of duality frames of our 3d Af = 2 theories (see [3] for physical explanation 
of this correspondence, in the case associated with 3-manifolds). This includes the case 
where 4d Af = 2 theory are complete in the classification of [60], and (in addition to several 
exceptional cases) the quivers coming from the triangulation, discussed in section 2.2. Our 

33 We can also consider loop operators located at the south poles, and the corresponding y-variables. This 
realizes the b ^ 1/6 symmetry of our theory. 
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analysis suggests that this correspondence is more general, and holds for non-complete 4d 
M = 2 theories; the classification program of the IR fixed points of 3d M = 2 theories is 
closely related with the classification of 4d M = 2 theories! 

5 Summary and Outlook 

In this paper, we have proposed a new conjecture that a class of 3d M = 2 theories are 
naturally and systematically associated with a sequence of quiver mutations. The data of 
the quiver mutation is encoded in a bipartite graph, the mutation network, from which we 
have computed the associated cluster partition function and identified the corresponding 3d 
M = 2 theory. The rules are summarized in Table 2. This is a rather general procedure, and 
includes in particular theories associated with 3-manifolds. 



Table 2: Dictionary between mutation network, 3-manifold and 3d M = 2 theory. 



mutation network 


3-manifold 


3d Af = 2 theory 


a white vertex 


an edge of tetrahedron 


an [/(l)-symmetry 


an intermediate white vertex 


an internal edge 


a gauge [/(l)-symmetry 


an initial/final white vertex 


a boundary edge 


a global [/(l)-symmetry 


a parameter associated 


a parameter on 


a vectormultiplet scalar 


with a white vertex 


an edge of a tetrahedron 




a black vertex 


an ideal tetrahedron 


an M = 2 hypermultiplet 


an edge connecting 


an edge belongs to 


U(l) charges 


black and white vertices 


an ideal tetrahedron 


of a hypermultiplet 


black vertices 


ideal tetrahedra 


a superpotential term 


connected to a white vertex 


glued around an edge 





We leave the detailed field theory analysis of our theories for future work. For example, it 
would be interesting to discuss the mutations of (A k , A n ) quivers discussed in [22, 61, 59, 62]. 

In this paper we have identified our 3d M = 2 theories based on the relations (1.1), 
(1.2) and the S% partition function. For the case with 3-manifolds, it is believed that our 3d 
M = 2 theories actually arise from compactification of 6d (2, 0) theory on a 3-manifold, and 
also from the boundary conditions of 4d M = 2 Abelian theories. The latter in particular 
gives a direct physical method to identify the 3d M = 2 theories, which are expected to have 
the same S% partition function as our 3d Abelian theories. This program has been carried 
out for 1/2 BPS boundary conditions in 4d M = 4 theory [30, 31], and more results in this 
direction will appear in the upcoming work [63]. 



3s 



There are also more mathematical questions to ask — our partition function defines a 
knot invariant, and it would be desirable to define the invariant more rigorously. In fact, the 
discussion in in section 3 uses braid groups and plat representations of knots, which is often 
used in the study of Jone polynomials and their generalizations (cf. [64]). 

For the case without 3-manifold description, we have different questions to ask. Why does 
the relation (1.2) hold? For precisely which class of 3d gauge theories does the relation hold? 
Do we have a string theory realization of our theories? As already commented previously, 
our discussion includes the case where our quiver is identified with the BPS quivers of 4d 
M = 2 theories, and for these examples it is expected that our 3d M = 2 theories are the 
1/2 BPS boundary theories of the 4d M = 2 theories. However, our quivers in this paper 
can be arbitrary quivers and they will not necessarily be the BPS quivers. 

The fact that (1.2) holds for a rather rich class of 3d Af = 2 theories is a strong indication 
that there exists a rather rich structure in the "space of 3d M = 2 theories" beyond the realm 
of 3-manifold theories, and what we know right now is probably only a tip of the iceberg of 
a much richer structure. Indeed, the cluster algebras our paper, and their interpretation as 
the algebra of loop operators, suggest that general philosophy that the IR fixed points of 3d 
M = 2 theories can be characterized by the algebra of 1/2-BPS loop operators. 

One important clue for this ambitious program is the mathematical structures discussed 
in this paper, such as cluster algebras and hyperbolic 3-manifolds. They appear in a diverse 
area of physics and mathematics, including wall crossing phenomena of 4d M = 2 theories 
[60, 65, 66], dimer integrable models [67, 68], on-shell scattering amplitudes [69], and super- 
potential conformal indices of 4d M = 1 theories and their dimensional reductions, as well 
as associated integrable spin lattices [70, 71]. 
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A Quantum Dilogarithms 

In this appendix we collect formulas for the so-called non-compact quantum dilogarithm 
function (simply called quantum dilogarithm function in the main text) s b (z) and e b (z) 
[27, 28, 29]. 

We define the function Sb(z) by 



s b (z) = exp 

and e b {z) by 



1 f 00 dw ( sin2zw z 



i Jo w \2smh(bw) smh(w/b) w 



(A.l) 



(\ f 00+10 dw e~ i2zw \ 

e b (z) = exp - / — , (A.2) 

V.4 > /_ 00+ io w smh(wb)smh[w/b) J 

where the integration contour is chosen above the pole w — 0. In both these expressions we 
require |Imz| < |ImQ,| for convergence at infinity. The two functions are related by 

e 6 (z) = e 2 e 24 Sfe ( z ) ; (A.3) 

where Q := b + b^ 1 . 

It immediately follows from the definition that 

e b {z) = e b -i(z) , s b (z) = s b -i(z) , s b (z)s b (-z) = 1 . (A.4) 

In the classical limit b — > 0, we have 

e b (z) exp (-2^Li 2 (-e 2 ^)) , (A.5) 

where 1^2(2) denotes classical dilogarithm function of Euler, defined by 

n=l " / ° 

The Fourier transform of e^z)* 1 basically comes back to itself [72]: 



dxe b (x) e 2mwx = e — 2 +Sd+Q 2 ) e h (w + i^ 



dxe b {x)~ l e 2mwx = e^-SKW) e& j _ w - z - 



QV 1 



(A.7) 
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B Hyperbolic Geometry in a Nutshell 



In this appendix we briefly recall minimal ingredients of classical 3d hyperbolic geometry 
(see, e.g., [53, 73]), for readers unfamiliar with the subject. Let H 3 be the 3d hyperbolic 
space, namely a upper half plane 

R 3 >0 = {(x h x 2 , y) | x x , x 2 e R, y > 0} , (B.l) 

with the metric 

ds 2 = ( dx i) 2 + ( dx i) 2 + d y 2 , B 2 ^ 

y 2 

An ideal tetrahedron is a tetrahedron whose all four vertices are on the boundary of H 3 
(see Figure 24). By a suitable isometry of H 3 we can take the vertices to be at positions 
0,1, z and infinity. This complex parameter z is called the modulus (shape parameter) of 
the tetrahedron. 








Figure 24: An ideal tetrahedron in H 3 has all the four vertices on the boundary of H 3 , which 
we can take to be {0, 1, z, oo} e C U {oo}. 

A tetrahedron has 6 edges and correspondingly 6 face angles. For an ideal tetrahedron 
with modulus z, the two face angles on the opposite side of the tetrahedron are the same. 
These are given as the arguments of three complex parameters 

z , z' := — , z" := l-z~\ (B.3) 

satisfying zz'z" = —1. These three distinct parametrization of tetrahedron modulus will 
play an important later, when we discuss our rules. 

When we glue the tetrahedra to construct 3-manifolds, we need to ensure that the sum 
of angles around an edge sum up to 2n. Since the angles of tetrahedra could be described 
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either by z,z' or z" depending on the parametrization, we have 

\\ *M Y[ z /(m) Y[ z"^ = 1 , (B.4) 

m: type 1 m: type 2 m: type 3 

where we classified a tetrahedron m depending on the angle around the edge is parametrized 
by z(' m \ z'^ or z"^ m K We call these equations gluing equations. 34 In the boundary torus 
this simply is represented the the condition that the product of z around a vertex is trivial. 
We will come back to this equation in section 3.3. 



C Dehn Twists as Flips 

As we explained in the main text, an element of the mapping class group could be represented 
as a sequence of flips on the 2d triangulation. In this appendix we identify this flip sequence 
systematically. This result will be important for practical computations. 

The mapping class group is generated by the Dehn twist _D 7 along non-contractible cycles 
(Figure 25). Moreover explicit generators and relations for the mapping class group of 
are known in the literature (see e.g. [52, 74]). This means that all we need to do is to identify 
an explicit sequence of flips corresponding to a single Dehn twist. 

G5B * (X/ 

Figure 25: A Dehn twist along a non-contractible cycle (colored red). 

Suppose that we perform a Dehn twist along a non-contractible cycle 7, which intersects 
several triangles as in Figure 26 (a). We can then verify that the flips shown in Figure 26 
(b) (c) realize the Dehn twist. Note also we need to exchange the labels appropriately after 
the flips. This is a rather general rule, which applies to any triangulation and to surfaces of 
any genus. 

The situation is especially simple in the case of the 1-punctured torus, studied in detail 
in [1, 13] — the Dehn twist along the a, /3-cycles both correspond to a single flip. 

For the n-punctured sphere, consider a cycle encircling the i-th and (i + l)-th punctures, 
and the corresponding Dehn twist Sj. These Dehn twists generate the braid group. We can 
verify that this Dehn twist Sj, in the triangulations in Figures 27 and 28, are given by either 
2 or 4 flips. 




34 They are also called structure equations. 
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Figure 26: A flip sequence for a Derm twist along a cycle 7. (a) The cycle 7 intersects many 
triangles along the way, some of which have an extra edge on the left and some on the right. 
We have highlighted all the triangles on left side (colored red) and some on the right (colored 
green), (b) For each red/green block, we perform a sequence of flips as in this Figure, (c) 
By repeating the (b), and we effectively shift the position of one triangle, and by repetition 
we obtain a Dehn twist on the triangulation. 
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Figure 28: The Dehn twist (sj) on the n-punctured sphere could be traded for two or four 
flips (and re-labeling). 
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